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This paper is a report of an updated measurement of the CP-violating phase tjig'^ and the decay width 
difference for the two mass eigenstates Ar s from flavor-tagged decay B® — ► J/ip<j>. The 68% confidence level 
intervals, including systematic uncertainties, are <f> J a ^^ = — Q.55^q'|| and AF S = Q.163^q'q|| ps — 1 . This 
measurement is in agreement with SM expected value, the p-value for the Standard Model point is 29.8%. The 
data sample corresponds to an integrated luminosity of 8.0 fb _1 accumulated with the DO detector using pp 
collisions at y/s = 1.96 TeV produced at the Fermilab Tevatron collider. 



1. Introduction 



In the standard model (SM), the light (L) and heavy (H) mass eigenstates of the mixed B® system are 
expected to have sizeable mass and decay width differences: AM S = Mh — Ml and Ar s = Tl — Th- The 
two mass eigenstates are expected to be almost pure CP eigenstates. The CP-violating phase that appears in 
^-i, b — > ccs decays, due to the interference of the decay with and without mixing, is predicted [l| to be <f> J J^ = 
-2(3 S = 2axg[-V tb V t * s /V cb V* s ] = -0.038 ± 0.002, where V tj are elements of the Cabibbo-Kobayashi-Maskawa 

~~ '> quark-mixing matrix New phenomena may alter the observed phase Q to <pi^^ = — 2j3 s + <j£. 

, The first direct constraint on <j) J s ^^ 0] was derived by analyzing B® — > J/ip(j) decays where the flavor (i.e., B® 

or B s ) at the time of production was not determined ("tagged"). It was followed by an improved analysis [Bj], 
based on 2.8 fb _1 of integrated luminosity, that included the information on the B® flavor at production. In 
-V ' that analysis we measured Ar s and the average lifetime of the P° system, T s = 1/T a , where T s = (T H + T L )/2. 
The CP-violating phase <\> J J^^ was also extracted for the first time. The measurement correlated two solutions 
for <f>i with two corresponding solutions for Ar s . Improved precision was obtained by refitting the results 
y—i ' using additional experimental constraints [6] . Here we present new results from the time-dependent amplitude 
t— I , analysis of the decay B° s — > J/ibd> using a data sample corresponding to an integrated luminosity of 8.0 fb _1 
collected with the DO detector [7| at the Fermilab Tevatron Collider. We measure Ar s ; the average lifetime of 

> ! the B° system, r s = 1/T S , where T s = (T H + T L )/2; and the CP- violating phase '. 

X: 



2. Data Sample and Event Reconstruction 



The analysis presented here is based on data accumulated between February 2002 and June 2010, corresponds 
to 8 fb _1 of integrated luminosity. 

We reconstruct the decay chain B° s — > J/ip<f>, J/ip — !• ~, — > K + K~ from candidate f J/ip,<fi) pairs 
consistent with coming from a common vertex and having an invariant mass in the range 5.37 ± 0.20 GeV. 
Events are collected with a mixture of single and dimuon triggers. To avoid a bias in the B® lifetime distribution 
events are rejected if they only satisfy triggers that impose a requirement on the track impact parameter with 
respect to the pp interaction vertex. 

P° candidate events are required to include two opposite-sign muons accompanied by two opposite-sign 
tracks. Both muons are required to be detected in the muon chambers inside the toroid magnet and at least 
one of the muons is required to be also detected outside the toroid. Invariant mass range for muon pairs is 
3.096 ± 0.350 GeV, consistent with J/tp decay. Jji\> candidates are combined with pairs of oppositely charged 
tracks (assigned the kaon mass) consistent with production at a common vertex, and with an invariant mass in 
the range 1.019 ± 0.030 GeV. Each of the four final-state tracks is required to have at least one SMT hit. 

A kinematic fit under the P|? decay hypothesis constrains the dimuon invariant mass to the world-average 
J ji\> mass 0] and constrains the four-track system to a common vertex. In events where multiple candidates 
satisfy these requirements, we select the candidate with the best decay vertex fit probability. 

The primary vertex (PV) is reconstructed using tracks that do not originate from the candidate P|? decay, and 
apply a constraint to the average beam-spot position in the transverse plane. We define the signed decay length 
of a B° s meson, L^ y , as the vector pointing from the PV to the decay vertex, projected on the B° s transverse 

momentum px- The proper decay time of a B° s candidate is given by t — Mb b L^ ■ p/(Px) where Mb s is the 
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world-average B® mass [8| , and p is the particle momentum. The distance in the beam direction between the 
PV and the B° s vertex is required to be less than 5 cm. Approximately 5 million events are accepted after the 
selection described in this section. 



3. Background Suppression 

The selection criteria are designed to optimize the measurement of <f>i^^ and Ar s . Most of the background 
is due to directly produced J /ip mesons accompanied by tracks arising from hadronization. This "prompt" 
background is distinguished from the "non-prompt" , or "inclusive B — > J/tp + X" background, where the J / ip 
meson is a product of a 6-hadron decay while the tracks forming the </> candidate emanate from a multi-body 
decay of a b hadron or from hadronization. Two different event selection approaches are used, one based on a 
multi-variate technique, and one based on simple limits on kinematic and event quality parameters. 

Three Monte Carlo (MC) samples are used to study background suppression: signal, prompt background, 
and non-prompt background. All three are generated with pythia [p. Hadronization is also done in pythia, 
but all hadrons carrying heavy flavors are passed on to EvtGen [lCj to model their decays. The prompt 
background MC sample consists of J ftp — > decays produced in gg — > J/ipg, gg — > J/tpj, and gj — ¥ J/ipg 

processes. The signal and non-prompt background samples are generated from primary bb pair production with 
all b hadrons being produced inclusively and the J ftp mesons forced into /x + yT decays. For the signal sample, 
events with a B® are selected, their decays to J/ip(p are implemented without mixing and with uniform angular 
distributions, and the B® mean lifetime is set to t s = 1.464 ps. There are approximately 10 6 events in each 
background and the signal MC samples. All events are passed through a full standard chain of GEANT-based [ll[ 
detector software of D0 simulation. 



4. Multivariate event selection 

To discriminate the signal from background events, we use the TMVA package [12]. In preliminary studies us- 
ing MC simulation, the Boosted Decision Tree (BDT) algorithm was found to demonstrate the best performance. 
Since prompt and non-prompt backgrounds have different kinematic behavior, we train two discriminants, one 
for each type of background. We use a set of 33 variables for the prompt background and 35 variables for the 
non-prompt background. 

To choose the best set of criteria for the two BDT discriminants, we start with 14 data samples with signal 
yeilds ranging from 4000 to 7000 events. For each sample we choose the pair of BDT cuts which gives the 
highest significance S/sqrt(S+B), where S (B) is the number of signal (background) events in the data sample. 
Figure []Ja) shows the number of signal events as a function of the total number of events for the 14 points. 
As the BDT criteria are loosened, the total number of events increases by a factor of ten, while the number of 
signal events increases by about 50%. 




Figure 1: (a) Number of — > J/ip<j) signal events as a function of the total number of events for the 14 criteria sets 
considered. (b) Mean value of (j(<p 3 ) as a function of the number of signal events. 



The choice of the final cut on the BDT output is based on an ensemble study. We perform a maximum- 
likelihood fit to the event distribution in the 2-dimensional (2D) space of B® candidate mass and proper time. 
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This 2D fit provides a parametrization of the background mass and proper time distribution. We then generate 
pseudo-experiments in the 5D space of B° s candidate mass, proper time, and three independent angles of decay 
products, using as input the parameters as obtained in a preliminary study, and the background from the 2D 
fit. We perform a 5D maximum likelihood fit on the ensembles and compare the distributions of the statistical 
uncertainties of <f> J s (er(</>s^*)) and Ar s (a(AT s )) for the different sets of criteria. The dependence of the 
mean values of o~{4> J s ^^) on the number of signal events is shown in FigureQJb). 

The mean statistical uncertainties of both <\> J J^^ and Ar s systematically decrease with increasing signal, 
favoring looser cuts. The gain in the parameter resolution is slower for the three loosest criteria, while the total 
number of events doubles from about 0.25 x 10 6 to 0.5 x 10 6 . The fits used for these ensemble tests were simplified, 
therefore the magnitude of the predicted uncertainty is expected to underestimate the final measured precision. 
However, the general trends should be valid. Based on these results, we choose the sample that contains about 
6500 signal events. 

We select a second event sample by applying criteria on event quality and kinematic quantities. We use the 
consistency of the results obtained for the BDT and for this sample as a measure of systematic effects related 
to imperfect modeling of the detector acceptance and of the selection requirements. The criteria are the same 
as in Refs. Q and @. We refer to this second sample as the "Square-cuts" sample. 



5. Flavor Tagging 

At the Tevatron, b quarks are mostly produced in bb pairs. The flavor of the initial state of the B® candidate 
is determined by exploiting properties of particles produced by the other b hadron ( "opposite-side tagging" , or 
OST). The OST-discriminating variables are based primarily on the presence of a muon or an electron from the 
semi-leptonic decay of the other b hadron produced in the pp interaction. If a charged lepton is not found, the 
algorithm attempts to reconstruct the decay vertex of the opposite-side b hadron and determine the net charge 
of particles forming the vertex. 

The OST algorithm, based on the Likelihood Ratio method, assigns to each event a value of the predicted 
tagging parameter d, in the range [—1,1], with d > tagged as an initial b quark and d < tagged as an initial 
b quark. Larger |d| values correspond to higher tagging confidence. In events where no tagging information is 
available d is set to zero. The efficiency e of the OST, defined as fraction of the number of candidates with 
d ^ 0, is 18%. The OST-discriminating variables and algorithm are described in detail in Ref. (l3j . 

The tagging dilution V is defined as V — A/" CO r — ■A/" wr /(7V CO r + A/" wr ) where N COI (N WI ) is the number of 
events with correctly (wrongly) identified initial B-meson flavor. 

The dependence of the tagging dilution on the tagging parameter d is calibrated with data for which the 
flavor (B or B) is known. The dilution calibration is based on four independent B® — > [ivD*^ data samples 

corresponding to different time periods. For each sample we perform an analysis of the B® — B° d oscillations 
described in Ref. [l4j. We divide the samples in five ranges of the tagging parameter \d\, and for each range we 
obtain a mean value of the dilution \T>\. The mixing frequency AMd is fitted simultaneously and is found to 
be stable and consistent with the world average value. The measured values of the tagging dilution \T>\ for the 
running period of time is parametrized by function: 

\V\ = ^ E2 (!) 

11 (l + exp(( Pl -M|)M)) (l+exp( Pl / P2 )) [> 

and the function is fitted to the data. There is a good agreement in the fits for running period of time, and 
hence a weighted average is taken. 



6. Maximum Likelihood Fit 

We perform a six-dimensional (6D) unbinned maximum likelihood fit to the proper decay time and its un- 
certainty, three decay angles characterizing the final state, and the mass of the B® candidate. We use events 
for which the invariant mass of the K + K~ pair is within the range 1.01 - 1.03 GeV. There are 104683 events 
in the BDT-based sample and 66455 events in the Square-cuts sample. We adopt the formulae and notation of 
Ref. [l53 ■ The normalized functional form of the differential decay rate includes an S- wave K K contribution in 
addition to the dominant "P-wave cj) — > K + K~ decay. To model the distributions of the signal and background 
we use the software library RooFit jl6| . 
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6.1. Signal model 

The angular distribution of the signal is expressed in the transvcrsity basis. In the coordinate system of the 
J/tp rest frame, where the </> meson moves in the x direction, the z axis is perpendicular to the decay plane of 
— > K + K~ , and p y {K + ) > 0. The transversity polar and azimuthal angles 9 and ip describe the direction of 
the positively-charged muon, while ip is the angle between p(K + ) and —p(J/ip) in the </> rest frame. 

In the transversity basis, the decay amplitude of the P° and B s mesons is decomposed into three independent 
components corresponding to linear polarization states of the vector mesons J /ip and <f>, which are polarized 
either longitudinally (0) or transversely to their direction of motion, and parallel (||) or perpendicular (_L) to 
each other. 

The time dependence of amplitudes Ai{t) and Ai{t) (i denotes one of {||,_L,0}), for B® and B s states to 
reach the final state J/ip (f> is: 

Ai(t) = F(t) [E+(t) ± e 2tf3 °E^(t)] a t , 

Mt) = F(t) [±E+(t)+e- 2tf3 °E_(t)] a t (2) 

where 

e -r s t/2 

F ® = 1 = , 0R( V (3) 

\JT H +T L ± COS 2(j s (T L - Tjj) 

and th and tl are the lifetimes of the heavy and light B® eigenstatcs. 

In the above equations the upper sign indicates a CP-even final state, the lower sign indicates a CP-odd final 
state, 

E±(t) = ^ ^(^^ +l ^) 4 ± e -\-^ +l ^> t \ , (4) 

and the amplitude parameters give the time-integrated decay rate to each of the polarization states, |oj| , 
satisfying: \cn\ 2 = 1. The normalized probability density functions Pb z (X z denotes one more same 

equation replaced by X z ) 

P Bz (0,y,i>,t) = Y^|A z (t) x n\ 2 , 

(5) 

where n is the muon momentum direction in the J/ip rest frame, 

n = (sin 9 cost/?, sin 9 sin ip, cos 9) , (6) 
and A(t) and A(t) are complex vector functions of time defined as 



m 

The values of Ai{t) at t = are denoted as Ai. They are related to the parameters a by 

l + (y-l)|aj 2 ' 

(8) 



\A,V = 



where y = (1 — z)/(l + z) and z = cos 2/3 s Ar s /(2r s ). By convention, the phase of A is set to zero and the 
phases of the other two amplitudes are denoted by 5\ | and <5j_ • 

For a given event, the decay rate is the sum of the functions Pb and P B weighted by the flavor tagging 
dilution factors (1 +T>)/2 and (1 — T>)/2, respectively. 

The contribution from the decay to J/tpK + K~ with the kaons in an S wave is expressed in terms of the 
S- wave fraction F$ and a phase S s . The squared sum of the V and S waves is integrated over the KK mass. 
For the V wave, we assume the non-relativistic Breit-Wigner model 
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g(M(KK)) = • M{KK) _\ u + ^ /2 (») 

with the (j> meson mass M<j> = 1.019 GeV and width I> = 4.26 MeV [8], and with AM(KK) = 0.02 GeV. 

For the 5-wave component, we assume a uniform distribution in the range 1.01 < M(KK) < 1.03 GeV. In the 
case of the BDT selection, it is modified by a iCfT-mass dependent factor corresponding to the BDT selection 
efficiency. We constrain the oscillation frequency to AAf s = 17.77±0.12 ps _1 , as measured in Ref. (l7| . TableQ] 
lists all physics parameters used in the fit. 



Parameter Definition 

|^4o| 2 P-wave longitudinal amplitude squared, at t = 

Al |A||| 2 /(1-|A,| 2 ) 
"f a (ps) B° mean lifetime 

AF S (ps -1 ) Heavy- light decay width difference 

F s K + K~ 5-wave fraction 

/3 S CP-violating phase ( = -<^ /V "V 2 ) 

h arg(A||/ylo) 

6 X a,Tg(A±/A ) 

$s arg(A s /A () ) 



Table I: Definition of nine real measurables for the decay B° — > J/ip4> used in the Maximum Likelihood fitting. 



For the signal mass distribution we use a Gaussian function with a free mean value, width, and normalization. 
The function describing the signal rate in the 6D space is invariant under the combined transformation /3 S — > 
7r/2 — j3 s , Ar s — » — Ar s , <5|| — > 2tt— 5i\, 5± — > tt — 8±, and S s — > n — S s . In addition, with a limited flavor-tagging 
power, there is an approximate symmetry around f3 s — for a given sign of Ar s . 

We correct the signal decay rate by a detector acceptance factor e(t/j, 9, ip) parametrized by coefficients of 
expansion in Legendre polynomials Pk(4>) and real harmonics Yi m (9, ip). The coefficients are obtained from 
Monte Carlo simulation. 



6.2. Background model 

The proper decay time distribution of the background is described by a sum of a prompt component, modeled 
as a Gaussian function centered at zero, and a non-prompt component. The non-prompt component is modeled 
as a superposition of one exponential decay for t < and two exponential decays for t > 0, with free slopes and 
normalizations. The lifetime resolution is modeled by an exponential convoluted with a Gaussian function, with 
two separate parameters for prompt and non-prompt background. To allow for the possibility of the lifetime 
uncertainty to be systematically underestimated, we introduce a free scale factor. 

The mass distributions of the two components of background are parametrized by low-order polynomials: a 
linear function for the prompt background and a quadratic function for the non-prompt background. The angular 
distribution of background is parametrized by Legendre and real harmonics expansion coefficients. A separate 
set of expansion coefficients cf m and cf m , with k — or 2 and / = 0, 1, 2, is used for the prompt and non-prompt 
background. A preliminary fit is first performed with all 17 x 2 parameters allowed to vary. In subsequent fits 
those that converge at values within two standard deviations of zero are set to zero. Nine free parameters remain, 
five for non-prompt background: c°_ l5 c^Jq, c^ 2 , Cq , and c\ 2 , and four for prompt background: c°_i, c 2 q, c 22 , and 
c\_x- All background parameters described above are varied simultaneously with physics parameters. In total, 
there are 36 parameters used in the fit. In addition to the nine physics parameters defined in Table |H they are: 
signal yield, mean mass and width, non-prompt background contribution, six non-prompt background lifetime 
parameters, four background time resolution parameters, one time resolution scale factor, three background 
mass distribution parameters, and nine parameters describing background angular distributions. 
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6.3. Systematic uncertainties 

There are several possible sources of systematic uncertainty in the measurements. These uncertainties are 
estimated for: 

• Flavor tagging: The nominal calibration of the flavor tagging dilution is determined as a weighted 
average of four samples separated by the running period. As an alternative, we use two separate calibration 
parameters, for the same running period. We also alter the nominal parameters by their uncertainties. 

• Proper decay time resolution: Fit results can be affected by the uncertainty of the assumed proper 
decay time resolution function. To assess the effect, we have used two alternative parameterizations 
obtained by random sampling of the resolution function. 

• Detector acceptance: The effects of imperfect modeling of the detector acceptance and of the selection 
requirements are estimated by investigating the consistency of the fit results for the sample based on the 
BDT selection and on the Square-cuts selection. Although the overlap between the two samples is 70%, 
and some statistical differences are expected, we interpret the differences in the results as a measure of 
systematic effects. 

• M(KK) resolution: The limited M(KK) resolution may affect the results of the analysis, especially 
the phases and the S- wave fraction F$, through the dependence of the S — V interference term on the V- 
wave mass model. We repeat the fits using this altered 0(1020) propagator as a measure of the sensitivity 
to the M{KK) resolution. 

The differences between the best-fit values and the alternative fit values provide a measure of systematic 
effects. For the best estimate of the C.L. ranges for all the measured physics quantities, we conduct Markov 
Chain Monte Carlo (MCMC) technique described in the next section. 



7. Confidence intervals from MCMC studies 

The maximum likelihood fit provides the best values of all free parameters, including the signal observables 
and background model parameters, their statistical uncertainties and their full correlation matrix. 

In addition to the free parameters determined in the fit, the model depends on a number of external constants 
whose inherent uncertainties are not taken into account in a given fit. Ideally, effects of uncertainties of external 
constants, such as time resolution parameters, flavor tagging dilution calibration, or detector acceptance, should 
be included in the model by introducing the appropriate parametrized probability density functions and allowing 
the parameters to vary. Such a procedure of proper integrating over the external parameter space would greatly 
increase the number of free parameters and would be prohibitive. Therefore, as a trade-off, we apply a random 
sampling of external parameter values within their uncertainties, we perform the analysis for thus created 
"alternative universes" , and we average the results. To do the averaging in the multidimensional space, taking 
into account non-Gaussian parameter distributions and correlations, we use the MCMC technique. 

The MCMC technique uses the Metropolis-Hastings algorithm [l8[ to generate a random sample proportional 
to a given probability distribution. The algorithm generates a sequence of "states" , a Markov chain, in which 
each state depends only on the previous state. 

To generate a Markov chain for a given maximum likelihood fit result, we start from the best-fit point x. We 
randomly generate a point x' according to the multivariate normal distribution exp(— (x' — x) • £ • (x 1 — x)/2), 
where £ is the covariance matrix. The new point is accepted if C(x')/ C(x) > 1, otherwise it is accepted with 
the probability C(x') / C(x). The process is continued until a desired number of states is achieved. To avoid a 
bias due to the choice of the initial state, we discard the early states which may "remember" the initial state. 
Our studies show that the initial state is "forgotten" after approximately 50 steps. We discard the first 100 
states in each chain. 

While we do not use any external numerical constraints on the polarization amplitudes, we note that the best- 
fit values of their magnitudes and phases are consistent with those measured in the U (3)-flavor related decay 
B° —> J/tpK* @, up to the sign ambiguities. Ref. [l|| predicts that the phases of the polarization amplitudes 
in the two decay processes should agree within approximately 0.17 radians. For S±, our measurement gives 
equivalent solutions near it and near zero, with only the former being in agreement with the value of 2.91 ± 0.06 
measured for — > J/ipK* by B factories. Therefore, in the following we limit the range of 5± to cos^ < 0. 
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7.1. Results 

The fit assigns 5598 ± 113 (5050 ± 105) events to the signal for the BDT (Square-cuts) sample. A single fit 
does not provide meaningful point estimates and uncertainties for the four phase parameters. Their estimates 
are obtained using the MCMC technique. 

Figure [U shows 68%, 90% and 95% C.L. contours in the ((pi^ ,Ar s ) plane for the BDT-based and for the 
Square-cuts samples. The point estimates of physics parameters are obtained from one-dimensional projections. 
The minimal range containing 68% of the area of the probability density function defines the one standard 
deviation C.L. interval for each parameter, while the most probable value defines the central value. 

The one-dimensional estimates of physics parameters for the BDT-cuts and Square-cuts sample are shown in 
TableHl 



Parameter BDT-cut sample Square-Cut sample Final Result 

~^7^) i.426±S:Sii ^IHl 

Ar s (ps" 1 ) 0.1291°:°™ 0.179™ 0.163t°:°g 

-o.49l8;3 -o.56l8:l§ -o.sslg:! 

\A \ 2 0.55218 8" 0.565 ±0.017 0.558±8'8w 

L4,,! 2 0.219ig^|? 0.249l8:8£ 0.23118^30 

5|| 3.15 ±0.27 3.15 ±0.19 3.15 ±0.22 

cos(<5 ± -<5 s ) -0.06 ±0.24 -0.20l8;2? -O.lllg;^ 

F s 0.146 ±0.035 0.173 ±0.036 0.173 ±0.036 



Table II: The one-dimensional estimates of physics parameters for the BDT-cuts sample, Square-cuts sample and final 
result values with systematic error. 




tfT* (rad) ftJW (rad) 



Figure 2: Two-dimensional 68%, 90% and 95% C.L. contour for BDT and Square cuts selection. The standard model 
expectation is indicated as a point with an error. 



To obtain the final C.L. ranges for physics parameters, we combine all eight MCMC chains, effectively 
averaging the probability density functions of the results of the fits to the BDT- and Square-cuts samples. 
Figure 1 shows 68%, 90% and 95% C.L. contours in the (^ J/# ,A1 S ) plane. The p-value for the SM point HJ 
((^ /W ,Ar s ) = (-0.038, 0.087 ps" 1 ) is 29.8%. 



8. Summary and Discussion 

We have presented a time-dependent angular analysis of the decay process B® — > J/ip<f). We measure B® 
mixing parameters, average lifetime, and decay amplitudes. In addition, we measure the amplitudes and phases 
of the polarization amplitudes. We also measure the level of the K K 5-wave contamination in the mass range 
1.01 - 1.03 GeV, Fs- The final result values for the 68% C.L. intervals, including systematic uncertainties, 
with the oscillation frequency constrained to AM S = 17.77 ± 0.12 ps -1 , are shown in the last column of Table 
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Figure 3: Two-dimensional 68%, 90% and 95% C.L. contours including systematic uncertainties. The standard model 
expectation is indicated as a point with an error. 



HH The p-value for the SM point (0s, Ar s ) = (-0.038,0.087 ps" 1 ) is 29.8%. 
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